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Properties of fluctuations in two chains of Tomonaga-Luttinger liquids coupled by the interchain 
hopping have been studied by calculating retarded response functions X/?(<fe, q y \ w) for charge and 
0^ ' Xa(<lx, q y \ uj) for spin where q x and q y (= or 7r) denote the longitudinal and transverse wave vector, 

respectively, and lj is the frequency. We have found the notable fact that the repulsive intrachain 
interaction results in the clear enhancement of lmx& (q x , tu) and the suppression of lmxf(q x ,Tr', tu) 
O |. at low energies. This result indicates the importance of the dynamical effect by the spin fluctuation 

with q y — n and small uj, which has a possibility to give rise to the attractive interaction for the 
^0 ' electron pairing. 
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I. INTRODUCTION 



| Two-chain system of Tomonaga-Luttinger liquids coupled by the interchain hopping and/or the exchange interaction 
r*j ■ is a basic model which connects one-dimensional interacting electron systems with quasi-one- and two-dimensional 
systems [[!]. When the intrachain interaction between conduction electrons is repulsive, the ground state exhibits the 
superconducting one [|2|-|9| where the spin gap found in the half-filling case [gjljj |l_2| still survives. The result is of 
interest for understanding the experimental fact that superconducting states have been observed in the doped ladder 
! systems, Sro.4Ca13.6Cu24O41.4s, under the pressure fl3[| . Further, the phase transition of quasi-one-dimensional organic 
^ • conductors, (TMTSF^X, which takes place from spin density wave state to superconducting one with increasing 
C*~) I pressure fll4|l5 could be related with crossover in the two chains [^6| 



It has been known that the superconductivity in the two-chain system results from the repulsive intrachain in- 
teraction of the backward scattering between electrons with anti-parallel spin [f7||l6[ when the interchain hopping is 
relevant p"7| . The study of the spin and charge fluctuations in terms of the response functions clarifies the dynamics 
00 ' of the two-coupled chains, and is useful for understanding the origin of superconductivity not only for two-chain 
^\ , systems but also for quasi-one-dimensional conductors. The response functions have been examined for two-coupled 
"j*^ ■ chain systems with only the forward scattering ]l8| ] , but those in the presence of the backward scattering are not yet 
clear. In the present paper, by use of effective Hamiltonian which is based on the renormalization group analysis, 
we investigate retarded response functions for charge and spin, given by X^(Qx> Qy', w ) and Xa(Qx, Qy', where q x and 
q y (= or 7r) are the longitudinal and transverse wave vector, respectively, and uj is the frequency. We demonstrate 
that the spin fluctuation dominates the low-lying excitation with q y = it from the calculation of the w-dependence of 
both hnx%(q x ,n;w) and ImXailx, tt; w). 
, In section II, the Hamiltonian for the two-coupled chain is represented by the phase variable and is investigated in 
terms of the new Fermion fields. The results of the calculation of the retarded response functions, Xp{lxi 1y \ u ) an d 
Xa(qx,<ly',u), are shown in section III. Section IV is devoted to summary and discussion. 

II. MODEL 



We consider a model of two chains coupled by the interchain hopping t where each chain consists of conduction 
electrons with repulsive intrachain interactions of both the backward scattering and the forward scattering. We note 
that the single chain in the absence of t leads to Tomonaga-Luttinger liquids Jl~9[ ] . The present system can be expressed 
by the phase Hamiltonian which is based on the bosonization method. In case that the energy is less than t, the 

present Hamiltonian is obtained by neglecting the non-linear terms including the misfit parameter 2qo = At/vp which 
is originated from the separation of Fermi wave vector due to the interchain hopping |^ Jl^j20| . Thus our Hamiltonian 
is given as, 
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where a' ~ «f/£ and Vp is the Fermi velocity. The phase variables 8± (</>±) and 8± (<j>±) express the fluc- 
tuations of the total charge (spin) and the transverse charge (spin) where [8 + (x),8^(x')] = [8 + (x),8-(x')] = 
[4>j r {x) 1 4>-(x')] = [<fi + (x), <fi-(x')] — i7rsgn(a; — x 1 ). The parameter g\ (#2) is a matrix element of the backward 
(forward) scattering, whose conventional definition is given by gi — > gi/(2nvp) flSfj . The coefficients in Eq.(|l|) 
given as vq 



are 
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order of Fermi energy. 

In Eq. (Q) , the excitation spectrum of the total charge fluctuation is gapless and has the same velocity as that in the 
absence of the interchain hopping. On the other hand, the interchain hopping plays important roles in the other degree 
of freedoms since the Hamiltonian expressing these fluctuation includes complex non-linear terms. For 32 — ffi/2 > 0, 
the renormalization to the strong coupling has been obtained for the terms proportional to cos \[28- cos \/2<j)+ , 
cos \/2(f>+ cos \/2(f>+ and cos V2~4>+ cos y/29- |l6|] . Thus the quantities, cos\/26L, cos\/20+ and cos\/20+ have finite 
expectation values, and the gapful excitations appear in the total spin, the transverse charge and the transverse spin 
degrees of freedoms. In case of the strong coupling, the behavior of fluctuations with low energy can be well described 
by use of the effective Hamiltonian in which the non-linear term with finite expectation value can be treated as the 
averaged one Following such a procedure, the parts describing the gapful excitation in Eq. (jl|) is rewritten as 

H' = ^( dxl-id^+f+^d^.)' 
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where (• • •) expresses the thermal average. In Eq.(||), the terms proportional to (^cos y/26-J cos \/2</>_ and 
(cos ^2^+) cos \/2(f>- have been discarded because these terms gives rise to the finite expectation value of cos V^-, 
which is inconsistent with the results derived from the renormalization group analysis. The quantities, ^cosy/29 

cos-\/2</>+/ and (cos\/2</>+) are related to the gap A s , A c and A s as 
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and are determined self-consistently similar to the previous case g\ = Jig]. 

Now we examine Eq. (|^) by introducing the Fermionic representation p2fl which has been applied for rju = 1 . 
Although the case rj^ = 1 is obtained as a special case in the calculation of the renormalization equation |l6[ , it is 
considered that such a setting does not change qualitatively the behavior of the solutions. The effect of 77^ 7^ 1 is 
discussed in the last section. Thus Eq.(0) is expressed as 
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H' = v F j dx {^(-ia^i) ~ 4Hd X ^)} +Asjdx + ^h] 

+ v F J dx ^ipl(-id x ip3) - ipl(-id x ip4,)^ + A C / dx {i0 3 04 - iV'lV'j} 
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where "01 ("02), "03 ("04) an( i 05 (06) are the field operators of right going (left going) Fermions corresponding to the 
transverse spin, the transverse charge and the total spin degree of freedoms, respectively. In terms of 0±, 6± and 0±, 
field operators, ipj are defined by 

\f2-KCt' 

^ b+n = 1 p^((- 1 )" < ^+ + '^-)pi(-l)' 1 f (^5 + jV6)+i7 r (jVi + W 2 + W 3 +jV 4 ) f ( 7 ) 

where n — and 1, and TV, (i = 1 ~ 6) is a number operator of the i-th Fermion. The Hilbert space is taken so that 
the numbers, N± + N2, N3 + Ni and N5 + Ng, are even integers ( JV{ (i = 1 ~ 6) : eigenvalue of iV{ ) in deriving 
Eq.(^|). Note that phase variables, 0±, 9± and <fr± can be also expressed as, 

i ± {x) = -Y,^z { - a ' W2+{qx) {D,{- q )±D 2 {- q )} , (8) 
8±(x) =-J2 ^e^/^ {D 3 (-q) ± D 4 (-g)} , (9) 
0±(x) = - 2 ^e(-'l«l/ 2 +^) {D 6 (-?) ± D 6 (-q)} , (10) 

9#0 q 

where Dj{— q) = J ipjipje~ lqx dx . 

In Eq.(o), the excitation spectra of the total spin, the transverse charge and the transverse spin degrees of freedom 



are respectively calculated as E k , s = y £|_ s + A;?, E k , c = y Cfc + Ac and £k ;S = y ^ + A 2 , where £fc,s = w^fc and 
£fc = vpk. By use of Eq.(||), the self-consistent equations Eqs.([|) - (||) for A s , A c and A s are rewritten as 
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By noting that the gap equations lead to A s = — A c = A due to <?+ > 0, we use E k = y £ 2 + A 2 = Ek, s = -E*,c in the 
following. In Fig.Q, the numerical results of the gap equations are shown with a choice of g\ = 0.45 and ta/vp = 0.1 
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FIG. 1. Solutions of the gap equations of Eqs.(ll) — (13) as a function of g+ — g*, for g\ = 0.45 and ta/vF = 0.1. The solid 
line and the dotted one express A s /£ c = — A c /£ c and — A 3 /£ c , respectively, where £ c is the cut-off energy of the order of t. 

III. RESPONSE FUNCTIONS 

We calculate the response functions denned by (y = p and a ) 



x (T T {v(x, 1; r) + th"v(x, 2; r)} {1/(3/, 1; 0) + e^i/fx 7 , 2; 0)}) 



(14) 

icj n — >w-\-i8 



where p(x, i;r) and a(x, i; r) denote operators of the charge and spin densities at the i(=l,2)-th chain, re- 
spectively. In Eq.(E3), the longitudinal wave vector q x is much smaller than 2k p ( kp is the Fermi 
wave vector ) and q y (= 0, 7r) is the transverse wave vector. Therefore the operators, v(l) ± f(2) (= 
u(x,1;t) ± i/(x,2;t)), in Eq.(14) are given as p(l) + p(2) = d x 9 + /n, a(l) + <j{2) = d x (t> + /n, p(l) - 
p(2) = {1/(27^0} Ep^P/^" tr(l) - a{2) = 

{l/(27ra')} J2 p a ii ppae lpm ° x exp{—ipp(9 + + p#_)/\/2}exp{ — ippa(<p + + p<p-)/y/2} where the summation denotes 

p = ±, a = ± and p = ± Jig ], The charge response functions Xp (fe, 0; u) with q a = 0, which is evaluated straightfor- 
wardly from the total charge density given by d x + /ir, is the same as that in the absence of the interchain hopping. 
On the other hand, Xa(lx,0',u), X^ilx, t; lo) and ^(c^, 7r; w) are calculated from the Fermionic representation of 
Eq.(|). These response functions at absolute zero temperature are given as 

X*( g .,0;o,) = ^( * + ., + '-—} , (15) 
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FIG. 2. The q x /go-dependence of the real part of the response functions, Rex P (qx,0;0), Rex<y{qx, 0; 0), Rex P {qx, 7r; 0) and 
^■ e X<y(lx, ?r; 0) which are normalized by 2/ttvf- Here the parameters are chosen as |A|/2t = 0.124 and |A s |/2t = 0.076 which 
corresponds to gi — 0.45, g-2 = 0.5 and ta/vF = 0.1. 

For u = 0, the real parts of Eqs.@-@ can be calculated analytically. The quantities, Rex P ((j , a: , 0; 0) and 
ReXo-(qx7 tt; 0), are given by 2rje/ r KVg and 2/ttvf, respectively, while Rexcr(<fa, 0; 0) and ReXp(o-) (<7z) tt; 0) are calcu- 
lated as follows, 



Re X ^,0;0)^ 1 /A 2 ^ 1^ + vW+4Ag| 

^ I v^qxy/iv^qx) 2 + 4A 2 - y / (v 4> q x )' 2 + 4A 2 | I 



2 v- 1 A 2 16 + V^+4A 2 i 
Re Xp ( g:c ,7r;0) = ]T - In i V - (19) 

In Fig.||, we show the real parts of the response functions which are normalized by 2/irvp- For q y = 0, Rexpfe, 0; 0) as 
a function of q x remains constant. The quantity ReXa^Qz, 0; 0) is reduced around q x — owing to the spin gap where 
the limiting behavior of Rexa(qx, 0; 0) for \q x \ <C | A s |/fJ<^ is given by 2/(i:v ( p) x u|g 2 /(6A 2 ). For q y — n, the gaps of 

the transverse fluctuations, A s and A c , lead to the suppression of Rexp(<7a;, tt; 0) around q x = ±qo. The minimum 
of Rexp(fe, 7r; 0) near q x = ±qo comes from the separation of Fermi wave vector which is caused by the interchain 
hopping. The quantity, Rex P (q x , 7r; 0), for \q x ± q \ < \A\/vf is expressed as 2/(ttvf) x {1/2 - (A/4i) 2 ln(4i/A) 2 + 
v F (q x ± go) 2 /12A 2 }. The quantity Rex CT ((Ja; ; ?n 0) as a function of becomes fiat due to the fact that the effect of 
A c and A s compensate each other in Eq.(pjj). 

Now we examine the w-dependence of the imaginary parts of Eq. ( |l7| ) which represents the spectral weight for charge 
and spin fluctuations with q y = n. In FigJ^, the quantities Yvax^{qxi 7r; u>) and lmx&(q x , 7r; u>), which arc normalized 
by 2/ttvf, are shown with the fixed q x /qo = 1-25 ( Fig.||(a) ) and 1.0 ( FigJ|(b) ). These quantities diverge at two 



locations given by uj L ^ H )/{2t) = y {1 — (+)q x /qo} 2 + 4{A/(2i)} 2 , which are ascribed to the separation of the Fermi 
wave vector by the interchain hopping. These singularities are also found in the absence of the interaction, where the 
corresponding response function becomes Xp ? ( C r)fe> 7r ) a ') ~~ * X R (Qx, tt; ui) with 

1 V" J vqx + qa , vq x + q \ 

X (9x,T;w) = — > ^\ -. r rr + -. r— — - \ ■ (20) 

v F {vq x + q ) - w - id v F {vq x + q ) + w + id 
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FIG. 4. The total weights, S^ (ty) and Sf (a) as a function of q x /qo- The dotted lines, S L and S H are those in the absence of 
the interaction. The parameters are the same as those in Fig. 2. 

In order to examine such a notable fact, we evaluate the weights by dividing Eq.([l7]) into two parts, i.e., 
Im Xp( CT )(9z, 7r ;w) = lmXp(^(qx,Tr;oj) + lmXp("(q x ,ir;uj), where hnx p £)(<lx,w;u) ( lm Xp(" t; w) ) comes from 
the part with z/ = — (+), and leads to the singularity near ui = ojl (u>h). In Fig.H, we show the total weight of 
Im Xp('i)(?x,7r; w) and Imx^fe, tt; w) which are defined by 

S p(°) = ^4fJ Q dujIm Xp(a)(<lx,ir,u) , (21) 

pOO 

S p{<T) = ^4fj o dujIm X p("{q x ,n;uj) . (22) 
G 



FIG. 5. The pairing in the momentum space (k x ,k y ) is shown schematically by the enclosed-dotted curve. The black ( 
white ) circles show the Fermi points of the upper ( lower ) band. The Fermi points with k y = —ir are equivalent to those with 
k y = 7T. The sign + and — expresses the fact that the order parameters of the superconducting state is out of phase. 

From the imaginary part of the response functions with q y — n, we concluded that among the fluctuations with 
q y = it, the spin degree of freedom is dominant compared to the charge one. The excitation of spin degree of freedom 
with q y = 7r and low energy seems to be closely related to superconductivity of the two chains. In the present 
superconducting state, the pair of the electrons is formed between two chains with in phase. This corresponds to the 
intraband pairing with out of phase in the momentum space (k x ,k y ) where the Fermi points are given 

by (±(£)_f + 9o/2), 0) for the lower band and (±(fcp — go/2), 0) for the upper band, respectively. The superconducting 
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state is explained in Fig.|5| where the pairing of the electrons around the Fermi points are expressed by the enclosed- 
dashed curves. By noting that Xp(o-) (<Zx, tt', iw n ) = (— l/V) dwlnr)^ % (q x , tt; w)/(iu; n — w )i it is considered that the 
fluctuations connecting these two kinds of pairs, Xp(a)(<lx,^",^n) with q x ~ ±go act as the attractive force [ p3[ . Since 
the spin degree of freedom in such fluctuation is dominant compared to the charge fluctuation, it is possible that 
the spin fluctuation with q y — tt results in the superconductivity of the two chains. Also in the organic conductors, 
(TMTSF)2X, it is argued that the pair is formed between the chains from the theoretical analysis j24| of the NMR 
relaxation rate of (TMTSF)2C104 |2j|. Therefore the origin of the superconductivity of the organic conductor seems 
to be the spin fluctuation with q y = ir. 
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